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Abstract. Motivated by recent experiments, we extend the Joanny-deGennes calculation of the elasticity 
of a contact line to an arbitrary contact angle and an arbitrary plate inclination in presence of gravity. 
This requires a diagonalization of the elastic modes around the non-linear equilibrium profile, which is 
carried out exactly. We then make detailed predictions for the avalanche-size distribution at quasi-static 
depinning: we study how the universal (i.e. short-scale independent) rescaled size distribution and the ratio 
of moments of local to global avalanches depend on the precise form of the elastic kernel. 

PACS. 68.35.Rh Phase transitions and critical phenomena 



1 Introduction 

Recent experiments on contact-line dynamics near depin- 
ning [lj have allowed for a direct measurement of the 
renormalized disorder correlator, and of the avalanche-size 
statistics. It was concluded in [1] that the precise form of 
the contact-line elasticity may have an important impact 
on the distribution of avalanche sizes. Thus, a detailed 
prediction for the former is needed. Surprisingly, a review 
of the literature reveals that this problem appears to have 
been treated only neglecting gravity 2|M or in presence of 
gravity only in the case where the equilibrium configura- 
tion of the interface is flat and horizontal, as for a vertical 
wall and a contact angle 9 of 90° 011], see Figs. [I] and [2] 

1 [2] introduce the capillary length L c ad hoc as a cutoff for 
the contact-line profile in response to a <5-like force. 





Fig. 2. The coordinate system for a wall inclined by (p with 
respect to the vertical: the coordinate y, not shown, is per- 
pendicular to the plane of the figure. The gray shaded volume 
is not filled by the liquid, leading to the gravity subtraction 
(second term in second-to-last line of eq. pb.) 

In the latter case, the elastic energy as a function of the 
Fourier-transformed height profile h q takes the form: 



£[h] = 



1 



dq 
2V 



= + 



(1) 



Fig. 1. The coordinate system for a vertical wall, 
air/liquid interface becomes flat for x 3> L c . 



The 



with k the inverse capillary length, and up to terms of 
order ft 3 . In other cases, interpolation formulas have been 
proposed [5]. 

There is numerical evidence [B] that the contact-line 
energy (JT|) does not correctly describe the shape of the 
contact line depinning from a single defect. Surprisingly, 
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the form of e q as given by ([!]) seems not to have been ques- 
tioned, rather the discrepancies were attributed to higher- 
order terms in h. The latter are indeed present [7j 8,9, 
even for 9 = 90° and a vertical wall (ip = 0), but for small 
perturbations they are subdominant w.r.t. the dominant 
term ([I I. 

In the first part of the present article (section [3]) we 
show that for contact-angles leading to a non-flat profile 
as plotted in figure [T] or for an inclined wall in figure [2] the 
elastic energy (JlJ cannot be used, but must be replaced 
by the more general form 



3in(0) cos(ip) (r 2 - 1) [t(r + <) + !] sin 2 (6>) 



sin(0 + 



t(r 2 + 3rt + 3t 2 - 1) 



(2) 



Then we discuss how this form can be measured from the 
contact line profile in presence of a single defect. This 
section is pedagogical and self-contained and can be read 
with no prior knowledge in either wetting or disordered 
systems. 

In the second part of the paper (section [4]) we examine 
the consequences of the form ([2 ) of the elastic kernel e q for 
the avalanche statistics of a contact-line at depinning, i.e. 
a contact-line advancing on a disordered substrate, when 
the fluid level is adiabatically increased. We calculate the 
local (i.e. at a given point) as well as global avalanche-size 
distribution, including the scaling functions. We work at 
quasi-static depinning using an over-damped equation of 
motion. The point of the present work is to study the ef- 
fects arising from the precise form of the elastic kernel. We 
therefore assume that details of the dynamics as well as 
non-linearities can be neglected in comparison. We use the 
methods introduced in [TO] to compute the avalanche-size 
distribution from the Functional Renormalization Group 
(FRG) theory of pinned elastic systems (see also Refs. jTT] 
for a pedagogical introduction and [T2"irn3] for early numer- 
ical tests of this theory) . We recall and analyze the formula 
given in these references, hence this section can be read 
with no prior knowledge of disordered systems; however 
understanding its foundations requires the knowledge of 
the above mentioned literature. Then we derive a general 
formula for the size distribution for an arbitrary elastic 
kernel, not given priorly, and compute universal ratios of 
local-avalanche-size moments versus global ones. Finally, 
we apply our general result to the contact line with the 
elasticity (p|. 



2 Model 

Consider a fluid in a semi-infinite reservoir, bounded by a 
planar plate, in presence of gravity. The plate is inclined 
by an angle — ir/2 < ip < ir/2 with respect to the vertical, 
see figure [2] We consider the coordinate system X, y, Z 
where Z is along the vertical direction, and the equation 
of the plate is Z = — X cot tp. The contact line of the 
fluid is parameterized by (h(y),y) along the plate, hence 



it is at {X — —h(y)smtp.y,Z — h(y) cos (p) in our coor- 
dinate system. The fluid occupies the space Z < Z(X, y) 
and X > — Zta,np, hence also X > — h(y) sin cp, where 
Z(X,y) is the height of the fluid-air surface. We choose 
the coordinate Z so that the reservoir level is Z = hence 
Z(X -> oo,y) = 0. 

The energy E is the sum of the fluid-air interface en- 
ergy, proportional to the area, and the gravitational en- 
ergy. It is a functional of Z(X, y), with explicit dependence 
on h(y). Its full expression reads: 



£[Z, h] 



= / dy 



E[Z, h] 



dX 



i + [d x z] 2 + [d y z} 2 + -K 2 z' 



X> — h(y) sin cp 



l 



o K / dy 



dX X 2 cot 2 Lp 



— h(y) sin ip 



dy I dX 1 

x>o 



cos 9 dy h(y) 



7 is the surface tension, and 

1 



(3) 



(4) 



defines the capillary length L c . We call £ the reduced en- 
ergy. The first term is the area of the fluid surface; the 
second is the gravitational cost in potential energy for 
bringing a fluid clement from infinity at level Z = 0, filling 
the reservoir up to height Z(X, y). The second line takes 
into account that for ip > and h(y) > (resp. ip < and 
h(y) < 0) the volume element < Z < — X cot ip (resp. 
—X cot ip < Z < 0) is actually not filled by the liquid; this 
is the gray shaded region on figure [2] Similarly for ip > 
and h(y) < (resp. ip < and h(y) > 0) there is a corre- 
sponding missing volume element which must be added. 
The third line is the subtraction of the surface energy of 
the flat profile h(y) = (i.e. Z = and X > 0), necessary 
to make the problem well-defined. The last term comes 
from the difference between solid-fluid 7sf and solid-air 
7sa surface energies which defines the equilibrium contact 
angle, denoted 9, via the usual relation 7 cos 9 = 7sa~ 7sf- 
The definition of the contact line implies the additional 
boundary condition 



Z {-h(y) sin ip, y) = h(y) cos tp 



(5) 



The energy of a given configuration h(y) of the contact 
line is obtained as E[h] := E\Zfr,h\ where Zh(X,y) min- 
imizes £[Z, h] at fixed h{y) under the constraint ([5|. By 
translational invariance, the minimum-energy configura- 
tion of the contact-line itself, i.e. the minimum of E[h], is 
attained for a straight line h(y) = ho where ho denotes 
the equilibrium height. 

In the next Section we compute the elastic energy of 
the contact line, i.e. E[h] to second order in its deforma- 
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tions, i.e. E[h] - E[h ] = E el [h] + 0{h 3 ) with 

E e \[h] = i J €qh-gh q . (6) 

Here h{y) = J q h q e lqv and we denote f :— J For a 
uniform deformation h(y) = h it takes the form: 

E c i[h] = ^m 2 (h - h ) 2 L y , m 2 = e q=0 , (7) 

which defines what we call the mass to, i.e. m 2 is the 
curvature of the parabolic well in which the contact line 
sits because of gravity. We then calculate e q , already an- 
nounced in eq. ([2|. 

3 Contact-line elasticity 
3.1 Model in shifted coordinates 

We start by introducing a more convenient expression for 
the energy of the system. The constraint X > —h(y) simp 
in the domain of integration is tedious to handle, so we 
introduce the function z(x,y) as 



z(x, y) := Z(X = x- h(y) sin ip, y) 



(8) 



It satisfies the same boundary condition 2(00, 0) = 0. We 
have also defined 



My) : = Kv) - h o 



(9) 



i.e. z is still the height along the vertical axis, but we have 
shifted the X coordinate so that the integration domain 
is 



X > Xq 



-ho sin ip 



(10) 



£[z,h] 



dy / dx 

y J x>x a 



z{x,y) 2 



+ V 1 + [d x z{x, y)] 2 + [d y z{x, y)+h'{y) sin tp d x z(x 1 y)] 2 
,2 



cos 2 t ( 9sin( ( 9 / My) 3 

6 Jy 

+(siny3 — cos 9) / h(y) . 



(13) 



To derive this result, we have used the relation 

J dx^\ + [d x zY + [d y zY- J dx\ 

X>-h(y) sinip ^>0 

= h(y)sm<p+ J dX^l + [d x ZY + [d y Z] 2 - 1 , 

X> — h(y) sin tp 

which was then rewritten in terms of x. 

The new height function obeys the constraint 

z(x ,y) = h(y) cos cp , (14) 

i.e. it is specified on the edge x = x . 

3.2 Zero mode and calculation of the mass 

Let us first consider a uniform displacement of the line 
h(y) — h. We can then restrict to y-independent height 
functions z(x,y) — z(x), and the reduced energy takes 
the form 

r r . K 2 

Ly 1 £[z,h}= / dx y/1 + z'{x) 2 - 1 + — z{xf 



G 



h 3 cos 2 ipsin tp + /i(sin <p — cos 9) . (15) 



There are additional constraints at the boundary, 

z(x — xq) — h cos tp , (16) 

and z = at infinity, h is arbitrary and not necessarily 
the preferred value at equilibrium noted Hq. 

The linear variation of the reduced functional £[z,h] 
around an arbitrary configuration (z(x), h), by (5z(x), Sh) 
can be written upon integration by part as 



5£[z, h]/L y 

dx 8z(x) 



Using that the derivatives, evaluated at X — x—h(y) sin ip, 
satisfy 

d x Z{X,y) =d x z(x,y) (11) 
d y Z(X,y) = d y z(x,y) + h'(y)sin(p d x z(x,y) , (12) 

one finds that the energy is now a functional noted E[z, h] = 
"f£[z,h] of z(x, y) and h(y) with 



x>x 



k 2 z(x) — d x 



z'{x) 



s/l + z'{xY 



+ 5z(xq) cos(9(xq) + <p) — — cos 2 <p> sin tp h 2 6h 
+ Sh(sm tp — cos 9) . 

We have defined 



cos{9{x) + tp) 
sin(9(x) + tp) 



z'(x) 



sf\ + z'{x) 2 

1 

^l + z'{x) 2 ' 



(17) 

(18) 
(19) 



hence 6{xq) is the contact angle at the wall and 9(x) + tp 
is the local angle with respect to the vertical. Note that, 
because of k 2 > 0, the profile decays and the boundary 
term at infinity does not contribute. The constraint on the 
boundary implies the additional relation 



Sz(xq) = Sh cos ip 



(20) 
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Consider now the function z(x) = Zh(x) which minimizes 
the energy at fixed h, i.e. 6h = and Sz(x ) = 0. This 
leaves the variations in the bulk, given in the first line of 
(17), leading to the stationarity condition 



z/L c 



K 2 z(x) = 

It can be integrated once, 

„2 



z"(x) 



(1 + z'(x) 2 ) 3 / 2 



%z{xf = 1 - / 1 
2 ^/l + z'{xf 



(21) 



(22) 



where the integration constant was fixed by considering 
the limit of x — > oo. This yields at x = xq, 





1.5 






\ 1.0 




0.5 





x/L c 



Fig. 3. The height profile z(x), in units of the capillary length 
L c , for 6 — 30° and tp — 0. Different contact-angles < 9 < 
7r/2 are obtained by moving the graph left/right. The case 
7r/2 < 6 < 7r is obtained by the reflection z — > —z. 



-h 2 cos 2 <p = 1 — sm(9(xo) + ip) 



= 2 sin" 



V 2 



+ 7: - t/4 



(23) 



using that 1 — sin a; = 2sin 2 (j:/2 — 7r/4). The sign of the 
root hcos (p must be opposite to the sign of 9(x )+ip— n/2, 
and similarly when solving for z in eq. ( 22 I . This yields 



hcos tp 



1 7T 

-2L c sin ( -(9(x ) + p--) 



(24) 



Integrating once more, we obtain the height profile z = 
Zh(x) in the inverse form, 



x{z) 




The integration constant was chosen to satisfy the con- 
straint ( 16 ). 



The reduced energy of the uniform deformation h(y) = 
h is thus a simple function £{h) — £ [z^, h], whose deriva- 
tive, £'(h) is easy to obtain. Indeed we can use the general 
variational formula (JTTJ) around (^(ar), h) setting Sz( x) = 
5hdhZh(x) with dhZh(xo) — 1 from the constraint (16 1. 
The bulk contribution is zero, due to the "equation of mo- 
tion" (21 1 satisfied by Zh{x). The bou ndary contributions, 
i.e. the second and third lines in (171), give 



(26) 



£'{h)/L y = cos Lpcos(0(xo) + ip) 



- h cos (p sin ip + sin p — cos ( 



Using eq. ( 23 1 this simplifies to 

£'(h)/L y = cosfl(a;o) 



(27) 



where 6{xq) is an implicit function of h, using (24 1. We 
first note that this allows to recover the usual condition 
for equilibrium, £'(ho) = namely that the local contact 



angle 9(xq) = 9, the equilibrium contact angle. The equa- 
tion determining ho is thus (24 1, with 9(xq) = 9, 



T 7T 

h Q cos(p= -2L c sin ( -(9 + p- -) 



(28) 



Inverting eq. ( 28 ) one finds from ( 27 1 away from equilib- 
rium 

£ '( h ) U 2 L 1 M o~ 
= ZiKCOS W \ 1 n z h z cos^ ip 

L„ V 4 

1 — 2^ 2r2 cos2 ^ sm f ~ cos ^ ■ (^) 

The mass for the zero mode is thus, after various simpli- 
fications 

2 j£"(h n ) 7«\/2 cos v? sin 

to = = — , . (30) 

L y y/1 + sm(9 + ip) 

An interesting special case is that of a flat interface, 9 + 
p=\. Then 

to 2 = 7Ksin 2 (6») . (31) 



3.3 Elastic energy for arbitrary deformations 

We now go back to the general case for h(y). The reduced 
energy functional £[z,h] given in (13) explicitly depends 
on h(y) . In addition h(y) also enters in the boundary con- 
dition (14 1. It is important to distinguish these two de- 



pendences. 

We call Zh{x,y) the profile which minimizes the re- 
duced energy functional £[z, h] with fixed h(y), hence with 
respect to bulk variations 6z(x, y) only, 



S£ 



Sz(x,y)\ hul 



■[x h ,h]=0 



(32) 



i.e. a partial (functional) derivative w.r.t. z only. 

We now expand around the equilibrium solution as 
follows: 



h(y) = h + h(y) 
Zh{x,y) = z Q (x) + z(x,y) , 



(33) 
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where Zq{x) := Zh {x) is the equilibrium profile deter- 
mined in the previous section. Since z is of order O(h), 
to compute the elastic energy we need to expand the min- 
imum energy to second order in h and z. Expanding (13 1 



to second order in the explicit dependence on z and h we 
have 

£[z h , h] - £[zo, h ] = 6 {1) £ + S^£ + 0(h 3 ) (34) 
5^£ = S^£ + S^£ + S^£ . 
The first term is the linear variation, 

8E r 



8^£= 6 ^[z M 
oz 



z + -jt-[zq, ho] ■ h 
on 



(35) 



whereas 8^ 2,pS> £ are second order variations computed be- 
low. Let us start with S^£: 



z(x,y) 



y Jx>x 



k 2 z (x) - d x 



z' (x) 



8z(x ,y) cos(9(x ) + tp) 



■ cos 2 tp sin <ph,Q h(y) 



+ / h(y) (sin tp — cos 8) 



(36) 



The first line is the bulk variation which vanishes since 



zo(x) satisfies the stationarity condition (21). The rest is 
the sum of the boundary variation in z and the variation 
in h. Thanks to the exact constraint 



z(xo,y) = h(y) C0Sl P 



(37) 



one can check, using the results of the previous section, 
that this sum vanishes identically if ho is the equilibrium 
value. Hence we have 8^£ = 0. Since z also contains 
subdominant h 2 contributions, this is a quite useful ob- 
servation, as we do not need to worry about them in com- 
puting the elastic energy. They drop out by virtue of the 
equilibrium condition. 

Let us now study the three second variations. The first 
one is 

(38, 



1 



5 {2,l) g= Z~ z . — M . z 

2 ozoz 



K 2 z(x,y) 2 

y,x>x 

The second one is 

S 2 E 



[d x z(x,y)} 2 



[d y z{x,y)\' 



[1 + Z' (x) 2 } 3 / 2 [l + 4(x) 2 ]V2 



8 (2 ^£ = h-—^{z ,h } 
ohoz 



(39) 



z' (x) 



ly,x>x Q ^l + Z Q {xy 



d y z(x,y) , 



where we have used h'(y) = h'(y) to alleviate the notation. 
The third contribution is 



6 (2,3a) £ + 5 (2,3b) S 



(40) 



where we write separately the second derivative with re- 
spect to the explicit dependence on h(y) of the surface 
energy, namely 

z' {xf 



2.:i« 



C 1 ■ 2 

8 = sin ip 
2 * 



ti(y) 2 - 
y Jx>x v 1 + z 'o( x ) 2 



and of the gravitational one, 
1 



-—k cos <p sin ip h 



Kyf 



(41) 



(42) 



To compute the first two contributions one needs to spec- 
i fy t he properties oiz(x,y) which follow from its d efin ition 
(33 1. The function Zh(x,y) must obey equation (32), for 



any h. Hence we can expand this equation order by order 
in h. For h = it yields again the stationarity condition 



(21 1 for zq(x); to first order it yields 
[zq, h \ + h ■ 



S£ 



5z6z(x, y) | bulk 



ShSz(x,y)\ hn 



lk 



[z ,h ] = 



which, in explicit form, yields the equation obeyed by 
z{x,y) (to linear order in h which is all we need here): 



d 1 d 1 d 2 

X [I + Z' (x) 2 ] 3 / 2 X [l + 4(x) 2 ]V2 V 



= sin tp h"(y) 



z' (x) 2 ]i/ 2 



z(x,y) 
(43) 



It must be solved with the boundary condition ( 37 1 and 
z(x = oo, y) — 0. 

Remarkably, this complicated looking equation, which 
depends on the profile zq(x) known only in the implicit 
form (25), can be solved analytically. First of all, using 



the stationarity equation (21 1 satisfied by Zq(x), one notes 
that 



z(x,y) = zi(x,y) := - simp h(y)z' (x) 



(44) 



is a particular solution of ( 43 1 which vanishes at x = oo 
and takes the value 

%i(x ,y) = sin</jcot(0 + <p)h(y) (45) 

at the boundary. The full solution can thus be written as 

z(x,y) = zi(x,y) + z 2 (x,y) , (46) 



where Z2(x, y) satisfies the homogeneous equation, i.e. (43 
setting the r.h.s. to zero. The boundary condition 
implies 



14 



Z2(x ,y) 



-My) • 



(47) 



sm(9 + ip) 

Both Z\ and zi vanish at x = oo. 

To solve the homogeneous equation, we go to Fourier 
space in y-direction, and introduce a new variable in in- 
direction. We thus look for the solution in the form 

~M*,y) = — ^!—, I e^~h(q)F s (S(x)) 



sin(# + tp) j,. 
S(x) = sm(9(x) + tp) = 
q = q/n 



1 



(48) 

(49) 
(50) 
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with S(xq) = sm(0 + ip). Deriving eq. (49 1 w.r.t. x, and 



using pTj) and (22) to express the result in terms of S(x) 



and to use the equation of motion for z\ (inhomogeneous) 
and z 2 (homogeneous). 



we obtain the rule for changing the derivatives, 



_ V2(l-5)v/S+T 
o x = k d s ■ 

The resulting equation for the function Fg(S) reads 
(fS + l)F q (S) + (1 - S)(7S 2 + S- 4)Ft{S) 
-2(5 - l) 2 S(S + 1)F''(5) = , 



(51) 



(52) 



Inserting z = z\ + z 2 into ( 56 I and using the equality 



(57 I we get 



S (2,l) £ + 5 (2,2) £ = _ 



h(xo,y) 



Z2(x ,y) 
[l + z> (x)*]V 



;d x z 2 (x,y) 



where sm(9 + ip) < S < 1. The constraint (47 1 implies the 
boundary conditions F q (sm(8 + ip)) — 1 and F q (l) = 0. 
After some search, the general solution of eq. (52 1 which 
satisfies Fg(l) — is found to be 



2 y y [i + ^p/s 

Z2(x ,y) 



X—Xq 



[1 + 4(X) 2 ] 3 / 2 



d x zi(x,y) 



■ sin tp 



v,x>x V 1 + z a\ x ) 



(58) 



9r(S) 



g r (sin(8 + p)) 
(s) ^ (l + ^T 



x [2r 2 + 3V2rVS+ 1 + 35+1 
where we denoted 

r := yr 



q 2 



(53) 



(54) 



(55) 



We now discuss simplifications. Firstly, the last ter m i n 
eq. (58) exactly cancels <5£( 2 ' 36 ); this is shown using (441. 



Secondly, from (44 1, (21) and (14 1, we obtain 



d x zi(x,y) 
[1 + z' (x) 2 } 3 / 2 



= — sin p cos ip n 2 h(y) hp . (59) 



We now evaluate the second variation. Consider first 



the sum of ( 38 1 and ( 39 1 . Using the equation of motion 
(43) for z, the combination <5( 2,1 '£ + ^5^ 2 ' 2 ^S can be inte- 



This shows that the second line, half the third line and 
c5( 2 ' 3b )£ cancel; the remaining half of the third line gives 
the first term reported in eq. ( |60l) below. The second term 
comes from the first line of (|58[), using (48), so we get 
finally: 



grated by part. Therefore, we obtain for the combination 

[«5(2,l)£ + I ( 5(2,2)£] + ^(2,2) S 



1 . 



z(x ,y) 



[1 + ^(a;) 2 ] 3 / 2 



d x z(x, y) 



2 sin(0 + ip) 

\ sin 2 6 sm{6 + <p) [ h q h- q F 9 (S(x))d x F 5 (S(x)) 



(60) 



■ sin <p 



y,x>x 



h"(y) , Z ' o{x) = z(x,y) 



(56) 



To compute the second term we use rule (51 1, where at the 



The last term has been integrated by part w.r.t. y^\ To 
continue, we note the useful equality 



zi(x ,y) 

y [l + 4(*) 2 ] 3 / 2 

Z2(x ,y) 



[1 + z' (x) 2 } 3 / 2 



d x z 2 (x,y) 

d x zi(x,y) 



end S must be evaluated on the boundaryi? = sin(6* + p). 
To compute the first term we use the value (23 1 for ho- 
This yields our final result for the elastic energy, 

Eci[h] = i J e q h q h- q , (61) 



ith 



z' (x) 



+ sin<P I h"{y) 

ly,x>x y/l + z' {x) z 



z 2 (x,y) , (57) 
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sin(0) cos(y) 



-1) [t( r + *) + l] sin 2 (6») 



<(r 2 + 3r< + 3i 2 - 1) 



which is a consequence of the two different ways to inte- 
grate by part 

[d x z 2 (x,y)][d x zi(x,y)] 



sm(6 + (p) + 1 



(62) 



L 



y,x>x 



k zi(x,y)z 2 {x,y) + 
[d v zi{x,y)][d y z 2 {x,y)\ 



One finds that e q is a scaling function of q/n which re- 
produces formula ( 30 1 for the energy of a uniform mode 



[1 + z^(z) 2 ] 3 / 2 



[l + ^(a:) 2 ]V2 



2 Note that the integration by part of d y produces no sur- 
face term. This can be made rigorous considering a periodic 
modulation h(y). We thus restrict here to functions that can 
be written as sum of periodic modulations, or have compact 
support. 



e q= o = m as computed in the previous section, and which 
behaves as e q s» Krysin 2 9\q\ for large \q\. 

When p + 9 = tt/2, the equilibrium shape of the inter- 
face is flat. Thus the elastic energy is expected to simplify. 
Indeed, it becomes 



S- =sin 2 (6»)Wl + 



(63) 
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qL c 



Fig. 4. fi7y/e<j as a function of qL c , for 
10°, ... , 90° (from bottom to top). 

My If 




y/L c 



Fig. 5. Sample profiles h(y), in units of //(fi/y), plotted as 
a function of y/L c , for = 40°, <p — 0° (bottom) and = 
45° (top). The profiles are the response to a force inside the 
green box \y\ < 8 — 0.3L C . The profile has to be continued 
symmetrically to the left. 
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12 



= 40°, ip = 0°, € J^y) 




0.2 



0.5 



1.0 



1.5 



2.0 



2.5 



3.0 



qL c 



Fig. 6. Top: Blue points: Reconstruction of e q , in units of tvy, 



as a function of q/ 'k, using (661 from the points of the bottom 
curve of figure [5] i.e. a vertical plate geometry tp = 0. Solid 
red line: the analytical result. One sees a numerical problem 
appearing at large q, at q — n/S, where S is the box-size in 
figure [5] However this is already far in the linear asymptotic 
regime (dashed line). Bottom: The same plot for smaller q. One 



3.4 The contact-line profile as a means of measuring e q 

Suppose the contact line is in force-free equilibrium. Then 
pull on it with a force per unit length: 



F(y) 



28 



\y\ < 6) 



(64) 



This leads in linear response (i.e. for the quadratic in h 
clastic energy we are using) to the following profile 



%) - / 



dq cos(qy) sm(q8) 
7T €„ qd 



(65) 



The limit of a (5-like force is recovered in the limit of S — > 0, 
which eliminates the last factor of sm(q8)/(qS). However 
the latter makes the integral convergent at large q. We 
have plotted for 8 = 40° and 8 — 0.3 two solutions on 
figure[5| one for ip — 0° (bottom) and the other for tp = 45° 
(top). We note that e q can be reconstructed from the 
profile as follows 



e q = f 



sin(qS) 
~^8~ 



dy cos(qy)h(y) 



(66) 



A reconstruction of e q starting from the blue points on the 
top of figure [5] is given in figure [6] 



4 Avalanche-size distributions 
4.1 The model 

We now study the case of a disordered plate, which is 
immersed with velocity v into the liquid reservoir. This 
is the geometry of the experiment described in Ref. pQ. 
We use some of the notations defined there: x denotes 
the coordinate along the contact line (denoted y in the 
previous section) and u(x) the height of the contact line 
along the plate, in the frame of the plate. We reserve the 
notation h(x) to the contact-line height measured in the 
laboratory frame, with the relatiorj^j 



h(x,t) = u{x,t) — w 



w 



vt 



(67) 



3 In the experiment of Ref. pQ , there is an additional slowly 
varying offset. 
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Even though it is a dynamical problem, it is useful to 
introduce the energy 



H[u] = I 
Jo 



dx 



[u(x) - w] 2 + V(x,u(x)) + SE[u) . (68) 



Here m 2 is the mass of the zero mode q = 0, and 5£J[it] := 
\ I q l e i ~ e o\ u q u -q the remaining part of the elastic en- 
ergy at non-zero wave vector q ^ 0. For the contact-line, 



m 2 = eq and SE[u] are both given in (|62|). We neglect pos- 
sible non-linear elastic terms 



The function V(x, u) is 
a random potential, whose derivative d u V(x, u) is short- 
ranged correlated, modeling the disordered substrate. 

The contact line is pinned by the disorder, but also 
trapped in the quadratic well with curvature m 2 . Advanc- 
ing the well-position w by immersing the plate with ve- 
locity v leads to a motion of the contact line, which we 
now study in the quasi-st atic limit, i.e. the limit of small 
v. From the energy form (68 1 one can derive, upon vari- 



ous assumptions about the fluid, equations of motion (see 
e.g. [4,7J and references therein). Our assumption here is 
that in the quasi-static limit, they lead to the same statis- 
tics as the simplest over-damped model studied in Ref. 
[TTTj . at least to the level of approximation that we use 
here (i.e. lowest order, i.e. one-loop FRG). Furthermore, 
although quasi-static dynamics and pure statics are dif- 
ferent, they lead, in the same order of approximation, to 
identical rescaled avalanche-size distributions. Deviations 
are expected only at the next, i.e. two-l oop order. This 
justifies our studying of the energy form (68), and, as we 



will see, our point is that noticeable effects already arise 
from the precise form of the elastic kernel. 



4.2 Global statistics of avalanches 

4.2.1 Definitions 

As can be seen on figure M the experiment [T] shows, as 
predicted by the theory (1211X0] . that the motion of the 





170 255 

w (urn) 

Fig. 8. Data from pQ: Height of the contact line h(w) averaged 
over 2L C , as a function of the position w of the plate (system: 
iso/Si). The fast depinning events (upwards) are clearly visible. 
Between them, the contact line moves downwards at the plate 
velocity v (here 1 fim/s). The straight line is the reference level 
ho- 



contact-line proceeds by sudden jumps, i.e. avalanche mo- 
tion. For an avalanche occurring at a given position w of 
the center of the quadratic well, one defines u~(x) the 
contact-line position just before the avalanche and u+(x) 
its position just after. The size of the avalanche is defined 
as 



S := dx [u+(x)-u w (x)] = / dx [h+(x) - h w (x)] . 
Jo Jo 

(69) 

Hence one can measure the avalanche size S as the area 
swiped both in the u or in the h coordinate system. 

We now recall the main results of Refs. 12,10 and 
apply the general formulation to various cases, including 
the contact line for which we have computed the elasticity 
in section 3. The characteristic rescaled function Z(A) is 
defined as: 



Z(X) := ^ (e**/*» - l) 

_ (S 2 ) -A'(0+) 
m -~ 2 (S) ~ to 4 



(70) 
(71) 



All averages (...) are over the normalized probability den- 
sity P(S) of avalanches. By definition Z(X) = A + A 2 + 
C(A 3 ). The function A(w) is the renormalized correlator 
of the disorder defined and studied in P~2~llTT7j and mea- 
sured in Ref. pQ. It is m dependent: at large m it is equal 
to the bare disorder correlator, while as m is decreased it 
develops a linear cusp at u = 0, whose value, A'(0 + ), is 
related to the second moment of the size distribution as 



Fig. 7. Sketch of the experimental setup used in [T]. The size 
of the image in the inset is 1.5 mm. 



displayed above in equation (71). 

The scale S m is the large-scale cutoff for avalanche 
sizes, originating from the quadratic well which suppresses 
the largest avalanches. It is an important scale as it allows 
to define universal functions in the limit where it becomes 
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large i.e. S m 3> S m in, where S m in is the typical size of 
the smallest avalanches. In the variable s := S/S m the 
avalanche-size distribution becomes universal. Universal 
means independent of small-scale details, but not of the 
large-scale setting, e.g. it will depend on the precise form 
of the elastic kernel. Indeed, one of the predictions of the 
FRG theory is that if the exponent r satisfies 2 > t > 1 
which will turn out to be the case here, then the distri- 
bution of avalanche sizes for S 3> S'min takes the form as 
m -> 0, i.e. S m > S'min, 



iJm \ ^ in J 



(72) 



The function p(s) is universal in the above sense. For a 
fixed elastic kernel it depends only on the space dimension 
d. Note that the normalized probability P(S) depends on 
the cut-off S m in via the first moment (S) which cannot 
be predicted by the theory, hence is an input from exper- 
iment. It is important to stress that while the function 
p(s) is universal and convenient for data analysis, it is not 
a probability distribution and is not normalized to unity. 
Rather, it satisfies from its definition (72 1 and using (71 1 
the two normalization conditions 



(s) p = J ds s p(s) = 1 
<* 2 ) p = / d SS 2 p(s) = 2 



(73) 
(74) 



Here and below we use the notation (s) p to denote an 
integration over p(s) and distinguish it from a true expec- 
tation value over P(S), denoted (...). Note that Z(X) and 
p(s) are related by the Laplace transform 



Z(X) 



/"CO 

/ dsp(s) [e Xs - 1] 
Jo 



(75) 



Note finally that the limit m — > means that L c is large 
compared to the microscopic cutoff a along the line, which 
may be of nanometer scale, or for strong disorder the size 
of the defects. 



4.2.2 Results from one-loop FRG 

We now summarize the results obtained in |10j for a gen- 
eral form of the clastic kernel e k . One defines: 



£fc/eo , 



eo = m 



(76) 



The general result of [TU] is that Z(X) satisfies, in an ex- 
pansion in powers of the renormalized disorder A(w) and 
up to terms of 0(A 2 ), 



Z = X + Z 2 + aJ{Z) , 



J{z) := 



'■Ik 



~-2 



Z 

(h - 2z) 



a=-e e- 2 A"(0+) (77) 

Jk 

9 

z z z 

2 e k -2z e k el 



(78) 



Note that the parameter 

e:=d c -d , (79) 

the distance to the upper critical dimension, introduced 
here for later purpose, cancels and is thus immaterial in 
this formula. As such, this equation is formally exact for 
any m, any e k and any dimension, up to 0{A 2 ) terms. 
Furthermore, consider now an elastic kernel such that: 

e k ~^oo KhP (80) 

with elasticity range (3 < 2 (i.e. long-ranged for (3 < 2 
and short-ranged for (3 — 2) and K the elastic constant. 
Assume that it takes the form 



h = e(K x l (i klpL) , n< 



P 



(81) 



where the elastic scaling function e(p) is constructed as the 
unique dimensionlcss function of the dimcnsionlcss argu- 
ment p which satisfies e(0) = 1, and e(p) = p@ for p — > oo. 
Then, the result of [10 is that the formula ( 77 1 can be con- 
trolled, in the limit m — > as an expansion in powers of 
e = d c — d, where d c — 2/3 is the upper critical dimension. 
For this one uses that 



e- 2 = p^K-W \ e(p)- 2 
Jp 

= ^K- 2+ i^L + o(\) 



(82) 



in the limit e = d c - d -> 0, where C fj = 2(Air)~ d ^ 2 /T{0) 
is independent of the shape of the function e(p). Since the 
1-loop FRG flow for A has the schematic form —fid^A = 



k c fe 



A 2 = (e f, e k 2 )A 2 , it is convenient to in- 



troduce the rescaled disorder, 
A"(0+) := (e / e^ 2 )A"(0+) = K- 2 C fjf ^A"(0+) , (83) 

Jk 

where the second equality holds to leading order in e. 
This is precisely the parameter a defined above, i.e. a = 
—A"(0). The one-loop FRG flow for the rescaled disorder 
admits a fixed point, from which one obtains, as m — * 0, 
that a flows to a = -A"(Q) = — |e(l - Ci)- Here Ci is the 
0(e) correction to the roughness exponent of the elastic 
object. For the type of disorder relevant for the contact 
line (i.e. random-field disorder) one has £i = 1/3, hence 



2 

'9 £ 



(84) 



This value will be used from now on in ( 77 1 , which is 
thus valid up to terms of order C(e 2 ). Using the above 
definitions one finds that to lowest order in e, J(z) depends 
only on the dimensionless scaling function e(p), and not 
on the stiffness K: 



J(z) - CI 



(e(p) — 2z) 2 e(p) — 2z e{p) e(p) 2 

(85) 

valid for SR elasticity as well as LR elasticity. Of course, 
since we work to lowest order in e, the above integral 
should be computed in the critical dimension d = d c . 

Note that the mean-field singularity of this equation 
corresponds to p — 0, e(p = 0) = 1, i.e. z = 1/2. 
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4.3 From Z{\) to p(s) 

In Refs. [T^IHU] we examined various choices for the elastic 
kernel e^, computed Z(\) for each choice, and extracted 
the scaled avalanche distribution p{s). Here we show that 
p{s) can directly be written as a function of Ck, or e(p). 



4.3.1 General formalism 



The definition ( 75 1 implies 

dssp(s)e Xs = Z'{\) 



(86) 



While p(s) does not admit a Laplace transform, the func- 
tion sp(s) does. Laplace inversion then yields 



sp(s) 



1 

2iTT 

1 

2in 



dXe~ Xs 



Z'(X) 



-too 
ioo 



dz e - s (Z-Z*-aJ(Z)) 



(87) 



up to terms of order a 2 . At this stage there is a heuristic 
step, to go from the A to the Z contour. We assume that 
the contour Z — ix with x £ K is the correct one, as it is 
for the mean-field case (a = 0) around which we perturb. 
We will check this result on known cases below. Then one 
has 

S P (s)=— dxe ^sx~sx 2 +sa,J(ix) (gg) 

2tt J_ 00 

/CO 
2 
dx J(ix)e~ %sx ~ sx 
-DO 

to lowest order in a. We have introduced the scaled mean- 



(89) 



field avalanche-size distribution, i.e. ( 87 1 at a = 0, 

For the elastic manifold, it holds for d > d c . We want to 
compute the correction to sp(s) of order a, i.e. of order 
e — d c — d: 



2tt 



dx J(ix]e 



isx — sx 



dt 



dx He-^-^dl 



k J — c 



+ e -t(e k -2 ix ) dy _ e -te kdy _ 3te -t~e kd 2y . „ ^ 



-isx—sx -\-iyx i 



dt 



dx xe 



-e^t— isx—sx 



x [e 2ltx (i-tx) +3tx-i] . 

(90) 

The integrand behaves (before integration over x and k) 
as t 2 at small t as a result of the counterterms. It is useful 
to introduce the elastic generating function, 



C(t) 



1 



1 



C/3 J p 



-te(p) 



(91) 



in terms of which we get, integrating (90 1 over x: 
p(s) 



Pmf(s) 



1 



1 



dt C(t)X(t,s) 



X(t,s) = e 



4i 3 + s 2 (2 + t) -2st(3 + 2t) 



2s + 3t(s - 2) 



(92) 



(93) 



This is our final and most general formula, valid for any 
elasticity and to first order in a, i.e. in e. By performing 
the integral over s, one checks that it satisfies automati- 
cally the two normalization conditions for p(s) (to order 
a). 



4.3.2 Special cases 

We now check that this formula recovers previous known 
results. 



Standard local elasticity: For e k = Kk 2 +m 2 one has (3 = 
2 and e(p) = p 2 + 1. Using that J k e" tfe2 = (Air)- d / 2 t- d ' 2 
and computing the momentum integral in d — d c = 4, one 
finds 



C(t) 



1 

2f 2 ' 



(94) 



Inserting into ( 92 ) and performing the t integral yields 

^L = l+^ [(ln S + 7E )( S -6)+4s-8^rV^+4] , 
Pmf(s) 16 

(95) 

which recovers the result (169) of [10] (to lowest order in 
a). 



Elasticity of flat contact line and generalization: From 
(63 1 the elasticity of a flat contact line ip + 9 = tt/2 is 
£fe = 7 sin 2 Q\Jk 2 + k 2 . This gives: 

m 2 = fj, = K/k (96) 



/3=1 , J ft:^7sin 2 6l , 77 

e(p) = ^p 2 + l . (97) 

It is instructive to slightly generalize this, and study 

e(p) - ^p 2 + (1 - a) 2 + a (98) 

with < a < 1, which interpolates between the flat con- 
tact line for a = and e(p) = \p\ + 1 for a = 1, two cases 
studied in Ref. [10J . From the definition (91 1, computing 
the integral in d = d c — 2, we obtain 



C(*) = —[! + (!-<*)*] 



(99) 



Inserting into ( 92 ) and performing the t integral yields 

P(s) 
Pmf(s) 



1 + | { 16 - 12a - 6 7E + (1 - a)v^s 3/2 



+ [(3 - 2a)s - 6] Ins - 4(3 - a)^^/s 
+ S (2a(5- 7E ) + 3( 7 E-2)} . (100) 
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For a = 0, this is the same as (E14) of [10], and for a = 1 
has the same form as ( 95 1 upon a rescaling of a by a factor 
of 2, as noted in Ref. |10j . 

To obtain nicer forms for p(s), and more convenient 
for extrapolations in physical dimension, one needs to re- 
exponentiate these direct e expansion results, as done in 
Ref. [TU] . We will not attempt to do this here for the gen- 
eral case. 



4.4 More on the generating function Z{\) 

Z(X) is easier to compare with numerical and experimental 
data than the disorder distribution p(s), since even if their 
statistics is mediocre averages over all data are taken, thus 
the statistical fluctuations are less pronounced. Therefore, 
we come back to the general formula and some examples. 



4.4.1 Generating function in the general case 

Since we work to the first order in a, i.e. in e = d c 
can write: 



d we 



~, s 1/ / \ J(|[l- V1-4A1) 



z—Z-k, 



J MF(A) 

(101) 

up to higher order terms, hence we only need to compute 
the function J(z) in eq. (78) at the point of the mean field 
solution z = Zmf(A) = |~(l - \/l - 4A) . 

The integral J{z) can be rewritten using the clastic 
generating function C(t) defined in ( |91[ ) for an arbitrary 
elastic kernel as 

/>oo 

J{z)= I dtC(t) [e 2tz (tz 2 + z) - z -3tz 2 ] . (102) 
Jo ' 



4.4.2 Standard elasticity 

Standard elasticity has the form e(p) — p 2 + 1 and using 
(94) one finds: 



J(z) = ^z [2z + (1 - 3z) log(l - 2z)\ . (103) 

This result is in agreem ent w ith eq. (150) of [10]. The 
value of a to be used in ( 101 1 for extrapolation to d = 1 
(i.e. e = 3) is a — —2/3. 



4.4.3 Flat contact line and generalization 



From the form e{p) — \/p 2 + (1 — a) 2 + a and using (99 1 
one finds 



J{z) 



2z(a + z — 3az) 
1 - 2z 



+ (a - 3z) log(l - 2z) 



(104) 



This agrees with (E10) of [TUJ for a = and with ([103) up 
to a global fac tor o f 2, as predicted in [TO]. The value of a 
to be used in ( 101 1 for extrapolation to d = 1 (i.e. e = 1) 
is a = -2/9. 



4.4.4 Contact line for arbitrary angles 9 and ip 

Having tested our general formula on the known case s w e 
can now apply them to the case of the elastic kernel ( 62 ) . 



Let us specify the relevant parameters and functions. For 
the contact line [3 = 1, d c — 2, and 



m 2 


= n = e q 


=0 — 


K 


= 7 sin 2 ( 






1 


[P9 


K 


" T D ~ 


V 7 


e g 


= e(p) = 


e(r) 


r 


= sji + 


*L = 

K 2 



7K 



■ sin u cos ip 



sm(6 + (p) + 1 



(r 2 -l) [t(r + t) + 1] sin 6 
(r 2 + 3rt + 3t 2 -1) cos if 

cos 2 (p K 
t 2 snr 9 n 



-Q ■ 



(105) 
(106) 

(107) 
(108) 
(109) 



We have defined for convenience a new function e(r) since 
r is the implicit variable in the final formula. Note that 
for the flat interface 9 + (p = tt/2, one recovers e(r) = r. 

The generating function Z(X) for the contact line can 
thus be computed from (101 1 using (85) and performing 
the integral in dimension d c = 2. Upon a change of vari- 
able we get: 



J(z) 



t 2 sin 2 9 

COS 2 if 



X I dr r 
l 



(110) 



: z~ 



[e(r)— 2z] 2 e(r)—2z e(r) e(r) 2 



Performing this integral analytically, or even the one in- 
volved in computing C(t), is rather awkward, and we have 
preferred to use numerical integration. 

To relate to the experiment of pQ, we have computed 
J(z) numerically for 9 = 40°, tp — 0°. We give some ex- 
plicit values: J(-l/2) = -0.493277, J(-l) = -3.10539, 
J(-3/2) = -8.79405, J(-2) = -18.1329. From this we 
have computed and plotted Z(X) in figure ^1 using the 
extrapolation to e = 1 in the one-loop result, i.e. using 
(101) with a = —2/9. This is compared to the result of 



the scaled kernel e(p) = \J p 2 + 1. The experimental data 
are also plotted on figure [9] Note that using the correct 
elasticity allows to get closer to the experimental data 
than using e(p) = \Jp 2 + 1, which is the kernel for the 
flat interface, and the only one previously available in the 
literature. 



4.5 Local statistics of avalanches 

4.5.1 Basic definitions 

It is also possible to study in experiments the statistics 
of the avalanches occurring within a given portion of the 
elastic object, e.g. of the contact line. One thus defines 
the size of a local avalanche, weighted by a characteristic 
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1.0 h 



0.5 




Fig. 9. The generating function Z(X) of avalanche sizes as de- 
fined in ( |75[ ). The curves are from bottom to top: Red/solid: 
mean-field result. Black dots: The experimental data of [T]. 
Dashed/green: Eq. ( |110|| int egrated numerically with the cor- 
rect elasticity kernel (l09\ for 6 = 40°, ip = 0°. Dot- 
Dashed/orange: The same with the elastic kernel for 6 = 90°, 



function 4>(x), e.g. a (5-distribution localized at x = 0, as 



4 : 



Sf>:= / d d x<t>(x)[u+(x)-u-(x)] 



& d x<j>{x)[h+{x)-h-{x)\ . (Ill) 



In |10j . we have shown how the characteristic function 



1 



AS* 



1 



(112) 



as well as the distribution of sizes can be computed within 
mean-field theory, valid at and above the upper critical di- 
mension d c . Since it already required some quite involved 
instanton calculus, we did not attempt to perform the e ex- 
pansion. Here we will perform the e expansion, to one loop 
0(e), not on the full distribution, but on a simpler quan- 
tity, the second moment. More precisely, we will compute 
the universal amplitude ratio between the second global 
and local moments, quantities easily accessible in experi- 
ments and numerics. In fact, it was already measured in 
[13] in the case of standard local elasticity and the present 
analytical result was quoted. Here we provide the details 
of the calculation and predict this ratio for the contact 
line which has not yet been measured. 

4.5.2 Universal amplitude ratios between global and local 
avalanches: general formulation 



From Eq. (F12) and (196) in [10] one has the general exact 
relation: 



— 



1 



xyzz' 



4>{y)(l>{x)g xz g vz , A zz ,[w] 



=0 + 



where g xy = J^^e," 1 and A zz \w\ = -R" zz ,[w\ is 
the functional derivative taken at a uniform w t = w of 
the functional R[w] (see e.g. [TUfTTj] for its definition), the 
limit w = + being taken at the end. For <p(x) — 1, using 
J zz , A ZZ '[w] = L d A'(w) one recovers 



(114) 



Here we use the notations of section (4.2.21 so that m 2 = 



We can now use the results of the e-expansion [15 ] [T6 [ 
[TTj . or its first-principle derivation in formula (463) of [TB] : 

A zz ,[w] = 8 zz ,A{w) ~ (g 2 zz , - 5 ZZ , Jf t )A>{w) 2 + 0(e 3 ) . 



Hence one obtains the ratio 

s <f, 

._ m 

a<t> — 



(115) 



f x H x ) Jxyz 

-2m 4 A"(0+) 



<t>{y) ( t ) i.x)g xz gy Z 
1 



<!>{y)<t>{x)g xz g xz >(gl z , - 5 Z 



xyzz' 



9t) + 0(e 2 ) , 



a formula valid for any function (f>(x). 

We know study avalanches which occur on a given sub- 
space of co-dimension dl ', defined by x_\_ = 0, i.e. of dimen- 
sion dcj, = d — d' , of the d-dimcnsional manifold, i.e. we 
choose 

<f>(x) = r)5 d '(x±) . (116) 
rj will be specified later. The above formula leads to 

= m A -qI 2 (d') 

+2mS [h{d)h{d') - I A (d, d')} A"(0) + ... 



2A"(0) 



I A (d,d') 
h{d)h{d>) 



(117) 



using the definition (83) of the rescaled disorder, A"(0) = 
A" (0)el2(d), and defining the integrals 



Analogous to ( 71 1 , we can define the scale of local avalanches 
as 

<(^) 2 ) 



:= 



2(5*) 



(113) 




(118) 



4 Note that the variable <l>(x) entering in 5"* and a? bears no 
relation to the plate angle ip introduced previously. 



d d 'k d d q 1 1 1 



(119) 
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Note that in the last integral the momentum in the upper 
bubble runs over a space of dimension d! < d. 

To perform the e expansion we now use the scaling 
form (81) of the elastic kernel. This yields to zeroth and 
first order in e = cL — d: 



= W d 'K- d '^I 2 (d') 



2A"(0) 



in terms of the dimensionless integrals 

d d q 1 




(2ir) d e(qy 



h(d,d') 
~h{d)~h{d') 

(120) 



(121) 



d d k d d q 1 1 1 



(122) 



(2Tr) d ' (27r) d e{k) 2 e(q) e(k + q) 
To obtain a meaningful amplitude we now make the choice 



V = 



(123) 



Note that comparing with experiments then requires a 
precise knowledge of the internal length scale K 1 ^ / fi de- 
fined by the (renormalized) mass m 2 = \iP and elastic 
constant K = lim 9S>K € q/l^ ■ Alternatively, it may be used 
as a method to measure this length. 



4.5.3 Standard local elasticity 



We know specify to standard elasticity = Kk 2 
i.e. p = 2, d c = 4 and e(p) = p 2 + 1. We study d' 
We need the integrals: 



= 1. 



/ 2 (d) = r(2- 

h{d' = 1 
I A (d,d<) 



(47T) 



1 



2 

dq 1 
27(<z 2 + l) 2 
I 



-d/2 
1 

4 



(124) 
(125) 
(126) 

2vr (k 2 + l) 2 J 2ir (2ir) d 

1 1 
[(k + q) 2 + p 2 + l](q 2 + p 2 + 1) ~ (q 2 +p 2 + iy 



h{d)h{d') h{d)h{d') 
This involves the integral 
j _ , dk 1 f dq d d - x p 



The calculation of this integral is performed in Appendix 
[A] To the order e we are working, we only need to compute 
it in d — d r — 4. We find 



i = (4tt)^5 r 



4+ -(7V3- 18)tt 



(128) 



Putting all these results together we obtain 

7tt ' 



O0 



1 

4 + 
0.25 



a 1 



6V3. 

0.0254934a + 0(e 2 ) 



0(e 2 ) 



(129) 



with a = — f (1 — Ci)- This is the expansion for a manifold 
with co-dimension one in e = 4 — d. 

For the 1-dimensional string, simulated in [13 , we need 
to extrapolate to e — 3 which yields a — —2/3. The two 
possible Pade expansions are 



1 
1 



1 



0.0169956 



= 3.74538 , 



16a 



1 



7 



1 7T 



.6^3 

Hence an estimate of the ratio is 

— = 3.74 ±0.01 



(130) 

3.72807 . (131) 
(132) 



Another expansion is the expansion in fixed dimension 
in powers of A. It is likely to be less precise than the 
e-expansion, but being simpler in spirit we indicate here 
its prediction. One evaluates directly the integrals in d = 
d' = 1. Using the results from appendix |A] this yields 



2 Z\"(0 H 
9 V 



(133) 



and inserting the one-loop value A"(0 + ) = —a = §e w e 



find the estimate l/a^ ss 3.81, slightly larger than (132) 



4.5.4 Contact-line elasticity 

For the contact-line elasticity, the e = 2 — d expansion 
can be done as indicated above. As it requires performing 
two-dimensional integrals we will only give here the fixed- 
dimension estimate for d = d' = 1, which requires only a 
one-dimensional integral on each momentum. We further 
perform them numerically using the form (109 1. We give 
two examples: 

(i) flat interface ip + 6 = vr/2. 

Then e(p) = \J p 2 + 1 and we obtain 



1 
2 

1 MM) 

/ 2 (1)/ 2 (1) 



0.312811. 



This yields 



1 + 0.312811 



(134) 
(135) 

(136) 



and w 1.87. 
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(ii) experiments of Ref. PQ, i.e. <p = 0°, 6 = 40°: 
We find 



J 2 (l) w 0.41867 



/ 2 (1)/ 2 (1) 



0.333101 ■ 



This yields 



0.41867 ( 1 + 0.333101- 



(137) 
(138) 

(139) 



and » 2.22. 

It is interesting to note that the "c lasti c length" r\ = 
K/ ft, which enters the definition ( 111 ), ( 116 1 of and 
I"G clcIs " 

k = f&*ir^+rt±± (140) 

cos y V 2 

i.e., it only involves the capillary length and geometric 
prefactors, thus is well-known experimentally. In addition, 
one can measure the dependence of on the plate angle 
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A Standard elasticity: Calculation of an integral 



The integral defined in the text can be written as 
dk f dq d d ~ x p 



I = 



2ttJ 2tv (2nY J tMM>0 



te 



((*+«) +P +l) + * 2 (9 +P +1) _ g-C*l+t2)C« +P + 



"] 



= (4tt)- 



*(tl+t 2 ) ^ [(«! + «2 + hh) V2 - (*(tl + t 2 )) 1/2 j 



1/2 | e -t-t!-t 2 



(141) 



/t,tl,t 2 >0 

using the general formula f a e~ q C q = (47r) _<i//2 (det C)~ 1//2 . Performing the changes of variables t\ — > tti and t 2 — > tt2 yields 

1 1 



(4»r)-^ f°dt /"dt! /°°df 2e - t(tl+t2+1 ^S-5(f 1 +^)5-| 
Jo Jo Jo 



V*2*i + U + 1 2 VU + te. 



Integration over t yields 



/ = (47r)-T-/ dti / dt 2 (*i +* 2 ) 3-3 (ti +t 2 + 1) 



JO JO 

Setting t\ — st, i 2 = (1 — s)t, one obtains 
I = 



2 2 



VMi + ii + i 2 Vti + 



(4^)-+ y o at y o ^-^t+ijT-r^-ij 

= (4?r) ~ y d« 1_ '(f + i)~r ( 



2 ~ 2 



Vl + (l-s)at 
2 



arccot I — = ] — 1 



One can perform the integral directly in d = 1, with the result 



(142) 



(143) 



(144) 



I = - 



144 



Thanks to the counterterm it also admits a limit for d = 4: 



/ = 



8tt 2 



I + £(-18 + 7^/5) 



Finally the expression for any d is 



V32 d r (2 - g ) 3*2 (§, 1, 2 - § ; -i §; |) + 20F (9V32 d - 43 d / 2 (d + 3)) r (f - |) 

4^+1^/29^3 



-2-d -d/2 



The combination 2 



[l - -'f^-'W 1 is plotted 

L i2(l)i2(d)J p 



on Fig. 
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(145) 



(146) 



(147) 
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